Equilibration of Luttinger liquid and conductance of quantum wires 
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Luttinger liquid theory describes one-dimensional electron systems in terms of non-interacting 
bosonic excitations. In this approximation thermal excitations are decoupled from the current 
flowing through a quantum wire, and the conductance is quantized. We show that relaxation 
processes not captured by the Luttinger liquid theory lead to equilibration of the excitations with 
the current and give rise to a temperature-dependent correction to the conductance. In long wires, 
the magnitude of the correction is expressed in terms of the velocities of bosonic excitations. In 
shorter wires it is controlled by the relaxation rate. 
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The conductance of quantum wires measured at low 
temperatures is quantized in units of 2e 2 /h [1]. This 
universality is commonly interpreted in the framework of 
the Luttinger liquid theory [2|, |3| , according to which the 
conductance is not affected by either the electron-electron 
interactions in the wire or the temperature p|. On 
the other hand, recent experiments show temperature- 
dependent corrections to the conductance , m appar- 
ent contradiction with the theory Q ■ While the origin of 
the corrections remains the subject of debate [lfj, much 
of the recent theoretical interest focused on the physics of 
one-dimensional fermions not captured by the Luttinger 
liquid theory }llMl5|. 

The concept of the Luttinger liquid applies to generic 
low-energy properties of systems of interacting one- 
dimensional fermions. The key feature of this theory is 
that the elementary excitations can be viewed as non- 
interacting bosons with an acoustic spectrum: 

H = J2 hv \lHh + ijW^ - ^o) 2 + vjJ 2 }, (1) 
1 

Ref . [2[ . Here b q is the operator destroying a boson with 
momentum Hq, the velocities v, vn, and vj may depend 
on the electron-electron interactions, L is the system size 
(periodic boundary conditions arc assumed). In addition 
to the boson occupation numbers N qi the state of the 
system is described by two integers, N and J. The latter 
can be related to the total numbers of the right- and 
left-moving electrons in the system, J = N R — N L , and 
N = N R + N L ; No is the total number of electrons in 
some reference state. 

It is important to note that within the Luttinger liq- 
uid approximation the numbers N R and N L of the right - 
and left-movers are conserved, i.e., electron-electron in- 
teractions do not cause backscattering of electrons. As 
a result, if the interacting quantum wire coupled to non- 
interacting leads is described by the Luttinger liquid 
model [3], the electrons entering the wire from the leads 
always pass through it, and the interactions do not af- 
fect the conductance. On the other hand, if the inter- 



actions are treated as a small perturbation to the non- 
interacting fermion picture, they do lead to backscatter- 
ing of electrons and reduce the conductance 13-111]. The 
backscattering occurs via three-electron scattering pro- 
cesses, which are not included in the Luttinger liquid 
theory. 

The corrections to the quantized conductance [H-Q are 
observed in the regime of low electron density, where the 
electron-electron interactions are strong. Thus the the- 
ory LJ, [15| based on the picture of weakly interacting 



electrons cannot be applied to the experiments directly, 
and the generalization to the case of arbitrary interaction 
strength is needed. In this paper we develop such a the- 
ory using the Luttinger liquid picture as a starting point 
and accounting for the electron backscattering processes. 

We start by noting that the Luttinger liquid Hamil- 
tonian ((T|) describes the stable low-energy fixed point of 
a one-dimensional system of interacting fermions. The 
full low-energy Hamiltonian consists of ([T]) and a number 
of perturbations describing interactions between the ex- 
citations. Although formally irrelevant, these terms are 
responsible for establishing the thermal equilibrium in 
the system. The equilibrium distribution of the bosonic 
excitations is dictated by conservation laws and indepen- 
dent of the detailed form of the irrelevant perturbations, 



e h(v\q\-uq)/T _ I ' 



(2) 



Due to the conservation of energy and momentum, the 
distribution function is characterized by two parameters, 
temperature T and velocity u. Relaxation towards the 
equilibrium distribution ([2]) is due to the irrelevant per- 
turbations which have the form of higher-order terms in 
bosonic variables, and the respective relaxation time is 
expected to scale as a power of temperature, tq oc T~ a . 

The boson distribution ^ can be derived from the 
general principle that for systems where energy and mo- 
mentum are conserved the probability of realization of a 
given microscopic state i depends on its energy Ei and 
momentum Pi as e ,-i E i— uP i)/ T . The momentum of a Lut- 
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tinger liquid is given by 

P = PF J + J2 h( lblb q , 



(3) 



where pf = ttHN/L and the first term represents the 
momentum of a filled Fermi surface with J = N R — N L 
extra electrons at the right Fermi point. The pres- 
ence of J in the momentum ^ means that its distri- 
bution e~( AJ ~ u pfJ)/t anc j -j-^e boson occupation num- 
bers ((2]) are controlled by the same parameter u. (Here 
A = Tthvj/2L.) In particular, the thermodynamic av- 
erage of J is given by J — Nu/vj. Taking into ac- 
count the relation Q between J and the electric current, 
/ = evj(J/L) 7 we conclude that at full equilibrium the 
velocity u of the boson gas coincides with the drift veloc- 
ity 



v d 



I L 
e N' 



(4) 



The conclusion u — Vd is obvious in the presence of 
Galilean invariance, where the boson gas should be at 
rest in a frame moving with the drift velocity. 

Establishment of full thermodynamic equilibrium re- 
quires the presence of backscattering processes which 
change the numbers of the right- and left-moving elec- 
trons N R and N L and transfer the corresponding mo- 
mentum to the boson excitations. Such processes are not 
accounted for in the Luttinger liquid theory and, as we 
discuss below, occur at a time scale r much longer than 
the time To required for the formation of the boson distri- 
bution (0). Therefore, at time scales t » r equilibration 
processes may be described by relaxation of the velocity 
u of the boson gas towards Vd , 



u-v d 



T 



(5) 



Due to the time dependence of u, the total momen- 
tum of the bosons changes as well. The rate of 
this change is easily obtained from Eqs. © and ([3]), 
Pb = (ttLT 2 /3Hv 3 )u. Conservation of the total mo- 
mentum (J3j) of the Luttinger liquid implies that the 
number of right-moving electrons changes at the rate 
N R = J/2 = -P b /2p F: resulting in 



T 2 



N*_ 

L 6H v 3 pf 



(6) 



Equations ([5]) and ([6]) describe equilibration of a system 
of interacting electrons at times t ^ tq in terms of the 
phenomenological parameter r. 

In the case of weak interactions the backscattering 
processes involve a hole diffusing through the bottom 
of the band due the three-particle collisions 15]. Be- 
cause of small occupation probabilities of the hole states, 
the backscattering rate follows the Arrhenius law, t _1 oc 



e -D/T ^ with activation energy D given by the Fermi 
energy. The momentum of the backscattered particle 
is redistributed among the particle-hole pairs near the 
Fermi points, which correspond to the bosons in the 
Luttinger liquid theory. In the case of strong interac- 
tions, the one-dimensional electrons form a Wigner crys- 
tal. The excitations of the crystal are phonons. At low 
energies they have linear dispersion and correspond to 
the bosons in the Luttinger liquid. The total (quasi-) 
momentum of the phonons can change via umklapp pro- 
cesses, whose rate also follows the Arrhenius law, with 
the activation energy D determined by the Debye fre- 
quency fiil |. 

The above examples demonstrate that the backscat- 
tering processes involve states at energies of the order 
of the bandwidth. Such processes are neglected in the 
renormalization group schemes describing the low-energy 
properties of Luttinger liquids. Although at arbitrary in- 
teraction strengths these processes have not been studied 
microscopically, we expect the scattering rate to follow 
the Arrhenius law with the activation energy of the or- 
der of the bandwidth, D ~ vp F - 

In the following we consider the conductance of a long 
uniform quantum wire adiabatically connected to ideal 
metal leads, starting with the case of spinless electrons. 
The interacting electrons in the wire form a Luttinger liq- 
uid, whereas the leads are assumed to be non-interacting. 
The numbers of the right- and left-moving electrons in 
the wire are controlled by the chemical potentials of the 
leads. In the case of a relatively short wire, L <C vt, 
the time of flight of an electron across the wire is shorter 
than r, and backscattering is negligible. In this case the 
conductance is Go = e 2 /h. In longer wires, L > vt, 
a significant fraction of right-movers are backscattered, 
and the electric current / through the wire is reduced, 



/ = G V + eN 1 



(7) 



Here V is the bias voltage. Equation ([7]) expresses the 
particle number conservation in the system. In the case 
of weakly interacting electrons it was derived in Ref . [3] ■ 
Since N R and N L retain their meanings as the numbers 
of right- and left-moving electrons in the Luttinger liq- 
uid, the derivation of Ref. 14[ applies at any interaction 
strength. 

We now use conservation of momentum and energy 
to find N R in Eq. ([7]) in the linear response regime 
I —± 0. Conservation of the momentum (j3|) implies that 
backscattering of a single right- mover, AN R = —1, is ac- 
companied by transfer of momentum 2pp to the bosons. 
On the other hand, according to Eq. (JIJ at J — > the 
total energy E of the bosons should remain unchanged. 
Thus, the momenta transferred to the right- and left- 
moving bosons must be equal, i.e., AP R L = pp and 
AE R,L = ±vpf- We therefore conclude that backscat- 
tering of right- moving electrons changes the energy of the 
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right-moving bosons at a rate 



E 1 



-vp F N' 



(8) 



In the absence of backscattering all the bosons in the 
system arc in thermal equilibrium with the leads, whose 
temperatures are assumed to be equal. Thus the energy 
current in the wire vanishes. Backscattering transfers 
energy from the left- to the right-moving bosons at the 
rate p]). resulting in the net energy current 



JE 



E 1 



(9) 



Equations ([T])-© express conservation of particle num- 
ber, momentum, and energy in the system. To make fur- 
ther progress we need to discuss the boson distribution 
function in a wire connected to leads. 

As argued above, an appreciable correction to the con- 
ductance arises when the wire is rather long, L > vt. The 
exponentially long relaxation time r describes establish- 
ment of the full equilibrium in the system. On the other 
hand, equilibration of the bosons to the distribution ([2]) 
occurs at the much faster time scale tq. Therefore at 
L > vt 3> vtq the boson distribution inside the wire 
takes the form (J2]) with velocity u that may or may not 
satisfy the full equilibrium condition u = Vd- In the lin- 
ear response regime J — ¥ only the bosons contribute to 
the energy current in the Luttinger liquid. Its evaluation 
for the distribution ^ is straightforward and gives 



JE 



ill 

3 hv 



(10) 



The parameters u and T of the boson distribution @ 
may depend on position x along the wire. In the linear 
regime, / — ¥ 0, one expects u oc / and dT/dx oc /. Then 
from conservation of energy, djp/dx = 0, and Eq. (|10j) 
we conclude that u is position independent. 

Combining Eqs. (|8l)- (fT0|) we find the following relation 
between the backscattering rate N R and velocity u, 



N 



r 2 



3 hv 2 pp 



(11) 



This equation expresses the energy and momentum con- 
servation laws in a wire connected to ideal leads and must 
be satisfied regardless of the equilibration rate r. On 
the other hand, if the equilibration rate is known, the 
backscattering rate N R can also be related to the boson 
gas velocity u via Eq. Comparing © and (fTl) we 
conclude that in a wire of length L the boson gas velocity 
u is related to the drift velocity (HI as 



Icq + L 



V d , ka = 2VT. 



(12) 



Substituting Eqs. CE2]) and <jTT]) into Q, we find the 
conductance of the wire 



G = G 1 



7T 2 T 2 L 



3 v 2 p 2 F l cq + L 



(13) 



The second term in Eq. (p~3|) represents a correction to 
the Luttinger liquid result [J] originating from equilibra- 
tion of bosonic excitations. The result (fT3|) is valid for 
any interaction strength. At weak interactions v coin- 
cides with the Fermi velocity vp = vp/ m, and Eq. ([13")) 
agrees with the results of Refs. 1 1 51 1 . In short wires, 
L -C l e qi the correction to the quantized conductance is 



proportional to the equilibration rate r 



- D ' T and 



Le- D < T . As the length 
the correction saturates 



the length of the wire, 5G oc 
of the wire increases, L 3> l e q, 
at 5G ~ (e 2 /h)(T/vpp) 2 , but remains small within the 
range of applicability of our theory, T <C vpp. Interest- 
ingly, the magnitude of the correction at L — > oo is fully 
characterized by the parameters of the Luttinger liquid 
Hamiltonian (TT]), with interaction strength entering via 
the speed of the bosonic excitations v. 

Our result (|T3|) applies to systems where electron spins 
are polarized by an external magnetic field. On the other 
hand, the temperature-dependent corrections to conduc- 
tance are observed in unpolarized wires, and spins 
are believed to play a significant role in this phenomenon. 
We now generalize our discussion to the case of electrons 
with spin. Compared to Eq. (fTJ, the Hamiltonian of 
the Luttinger liquid with spins [3j has twice the number 
of modes. Most importantly, there are two branches of 
bosonic excitations, corresponding to waves of charge and 
spin densities, and propagating at different velocities, v p 
and v ai respectively. In addition, instead of the numbers 
of right- and left-moving electrons N R ' L — (N ± J)/2, 
the state of each system is described by the numbers of 
such electrons for each spin direction, N^' L and N^' L . 

The discussion leading to Eq. (fT3")) can then be read- 
ily modified as follows. The conductance Go of non- 
interacting wire in Eq. ([7]) doubles in the presence of 



The number of right-movers N R in Eqs. ©, 
, and ([TT]) should be understood as Nl* + Nf 1 , 



spins. 
0, I 

whereas E R in Eqs. <j8j) and ([9]) becomes the total en- 
ergy of right moving bosons, in both charge and spin 
channels. In the right-hand sides of Eq. ([6]) and (fTQ|) 
one has to account for two types of boson excitations by 
replacing v~ n — > + v~ n . The most significant revi- 
sion appears in Eq. JSJ). When an electron is backscat- 
tered and momentum pp is transferred to right- moving 
bosons, it is distributed unevenly between the charge and 
spin branches. Due to the relatively fast equilibration of 
bosons their distribution functions have the form ^ with 
different velocities v p and v a in the two channels but the 
same u. Momentum change is accommodated by a small 
change of u. As a result, the fraction of the momentum 
pp transferred into charge or spin branch is proportional 



to v 



-3 

p,(7 5 



whereas the corresponding energies scale as v p 



We thus generalize Eq. (|5|) as [17| 



E R = 



-p F N 



R 



(14) 



Vp + V a 

As a result the equilibration length l cq in Eq. (jT2"j) be- 



4 



(15) 



and we recover the general form (IT3l) of the dependence of 
conductance on the length of the wire with the following 
modifications: 



2e 2 



(16) 



The presence of both v p and v a in the expression for the 
correction to the quantized conductance 2e 2 /h reflects 
the fact that the momentum of backscattered electrons 
is distributed among the bosonic excitations in both the 
charge and spin channels. 

The two velocities v p and v a coincide with the Fermi 
velocity vp in the limit of weak interactions. For repul- 
sive interactions, v p > vp and v„ < vp. The experimen- 
tally observed temperature-dependent corrections to the 
conductance of quantum wires are most prominent 
at low electron density. In this regime the Coulomb in- 
teractions between electrons are effectively very strong, 
e 2 /Hvp » 1, and v p ^> v a . The expression for the con- 
ductance of a quantum wire given by Eqs. (|13[) and (|16[) 
then simplifies, 



2e 2 



T 



L 



1- — 

6 V v a p F ) l cq + L 



(17) 



with Z cq = 2v a T, see Eq. (jl5[) . Importantly, the magni- 
tude of the correction to the quantized conductance is 
controlled by the spin velocity. The reason is that at 
v a <C v p both the energy and momentum of the bosonic 
excitations are dominated by the softer spin branch. 

Similarly to the case of spinless electrons, for long 
wires, L 3> l C q, the correction shows quadratic temper- 
ature dependence, SG ~ (e 2 / h)(T / Vo-Pf) 2 ■ Within the 
range of applicability of our theory, T <C v a pp, the cor- 
rection remains small. In short wires, L <C l eq , the tem- 
perature dependence of SG is controlled by that of the 
relaxation time, SG oc t . We are not aware of micro- 
scopic calculations of the relaxation time r for strongly 
interacting electrons with spin. However, based on the 
analogy with the case of spinless electrons [l6| , we expect 
activated temperature dependence r _1 oc e~ D °/ T , with 
the activation energy of the order of the bandwidth of 
the spin excitations, D a ~ v a pp -C v p pp. 

In a wire of fixed length the crossover between the 
regimes of long and short wires can be explored by chang- 
ing the temperature T. Because of the exponential tem- 
perature dependence of r, the condition L ~ 2v a T is 
satisfied at T ~ D a /\nN, where N is the number of 
electrons in the one-dimensional part of the device. At 
T < D a /\nN the activated behavior SG oc e^ D "/ T 
should be observed, in agreement with the experiment 



[6j. The quadratic temperature dependence is expected 
in the relatively narrow range D a j In AT <C T <C D a , and 
its experimental observation requires working with rather 



long wires. At the upper limit of this range SG 



At this point the system crosses over into the so-called 
spin-incoherent regime. Study of this crossover would be 
very interesting from both the theoretical point of view 
and because of its possible relevance for the experiments 
showing a shoulder-like feature in the conductance 
at 0.7 x 2e 2 /h. 

To summarize, we developed a theory of conduc- 
tance of quantum wires with strongly interacting elec- 
trons. Our theory extends the Luttinger liquid picture 
to account for the equilibration processes. Such pro- 
cesses involve backscattering of electrons and result in a 
temperature-dependent correction to quantized conduc- 
tance of the wire. 
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